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Abstract 

Solutions of the field equations of the Nonsymmetric Gravitational Theory with gr 0i n = are obtained for 
the homogeneous, plane-symmetric, time-dependent case, both in vacuum and in the presence of a perfect 
fluid. Cosmological consequences include a dependence of the speed of light on its polarisation, as in a 
birefringent crystal. 
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1. Introduction 



As an alternative to Einstein's theory of gravitation (GR), the Nonsymmetric Theory of Gravitation 
(NGT) possesses a richer structure with the potential for significantly different physical predictions. This 
structure, first proposed by Einstein and Straus (1945, 1946) in an attempt to unify gravitation and elec- 
tromagnetism, is characterised by a nonsymmetric fundamental tensor g. The tangent bundle of space- 
time is endowed with a hypercomplex algebraic structure (Crumeyrolle 1967, Kunstatter et al 1983, Mann 
1984, 1989, Moffat 1984) compatible with g and with the connections. In terms of components, one has 
Suv = + e £[iH' wnere e2 = +1' an d round and square brackets stand for symmetrisation and antisym- 
metrisation, respectively 

Although this attempt at unification was not successful, Moffat (1979) suggested that the extended 
structure could be interpreted as a generalised theory of gravitation. For several years thereafter, the 
physical content of the theory was thought (Moffat 1991) to arise from the existence of a new conserved 
charge which entered in the time-space components (g[ io j) of g, as well as in its diagonal ones (line element), 
but with the space-space components (g^j) being put equal to zero. After it was pointed out (Damour 
et al 1992, 1993) that with this choice, the theory, unprotected by gauge invariance, would allow coupled 
negative-energy modes to propagate to infinity, a version with g^-, ^ was proposed (Cornish and Moffat 

1994, Legare and Moffat 1995, Moffat 1995a, 1995b, 1995c) which avoids this defect, as well as other stability 
problems (Clayton 1996, 1997), as long as the gj io j can be put equal to zero. Claims (Moffat and Sokolov 

1995, Moffat 1995c) and counterclaims (Burko and Ori 1995) have been made since then about the absence 
of apparent event horizons under dynamical conditions in the spherically-symmetric case. Unfortunately, 
only the static vacuum spherically-symmetric solution has been obtained analytically (Wyman 1950) for 
&[M] = 0> an d if is difficult to infer from it the behaviour of time-dependent solutions. 

This paper presents the first such analytic dynamical solution, but instead of the spherically-symmetric 
case, we address the more tractable homogeneous plane-symmetric case which may have applications to 
cosmology. Section 2 sets the relevant field equations to be solved. Very little theoretical background is 
provided as this is easily accessible in the existing literature (see for instance Moffat 1995b). Section 3 
presents the solution to the vacuum case, and Section 4 the solution to the case with matter. Section 5 
contains a few comments about the matter solution, in particular expressions containing NGT corrections 
for the GR prediction for the rate of expansion of the universe. In Section 6 we derive one major consequence 
of the direct coupling of NGT to electromagnetic fields: the universe can behave somewhat like a birefringcnt 
crystal so that the polarisation vector of light is rotated as it travels over cosmic distances. The paper ends 
with a few concluding remarks in Section 7. 

2. Field Equations 

In this paper, boldface symbols denote tensors, while a tilde above a tensor makes it a density. We work 
in units where G = c = 1. Finally, X = X t = dX/dt. 

In the homogeneous, plane-symmetric case, we use the fundamental tensor: 

g = e v dt ®dt- e x dx <g> dx - (3{dy <g> dy + dz <g> dz) + f(dy A dz), (2.1) 

2 



where, without loss of generality, the yz plane has been chosen as the plane of symmetry, and v, A, (3 and / 
are functions of time. 

Denoting by g^ u the contravariant components of g, the homogeneous, plane- symmetric connection 
components obey 

?",x + s pv r; x + r p n P r v r[ Xp) = o (2.2) 

Define the Einstein tensor, G = R — \g TrR, where the covariant components of R are given by 

p pa _ 1/ra 1 pa \ pa p/3 _i_ pa /o a\ 

Then, the other homogeneous, plane-symmetric field equations reduce to 

G = 8ttT (2.4) 

where a possible contribution from a cosmological constant has been ignored. 

We assume matter to consist of a perfect fluid made of particles moving with four-velocity u. The 
energy-momentum tensor takes the form (Vincent 1985): 

T = (p + p)u®u-pg, (2.5) 

with p the internal energy density andp the pressure. We choose a comoving frame, where u l — (i — x, y, z), 
and u = 1. Since only three of the four equations (2.4) will be independent, one of the four functions in 
(2.1) is arbitrary; it makes sense to take v = 0. The time t is interpreted as the proper time of the particles 
in the fluid. 

Now, introduce the parametrisation (3 — R 2 cos?/;, and / = R 2 simp, as well as the definition A = \ogR 2 . 
The nonvanishing connection components which solve (2.2) read 



r* = 

XX 


±Ae A 

2 5 


yy 
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1 \ 

2 ' 


1 [yz] 


1 (tz) - 




r v - 



-q tv] = \i> (2-6) 

With the help of (2.3) and (2.6), (2.4) can now be written down in terms of the components of g. Indeed, 
the nonzero mixed components of G are 

G* =\{A 2 -ii 2 ) + \A\ = ^p (2.7) 

G x x =A+\(3A 2 +iP 2 ) = -8irp (2.8) 

c y y = g; 

= LA+\(A 2 +iP 2 ) + \A\+\\ 2 + ±\ 

= -8irp (2.9) 

G y z = ?A + i^ + iA^ = (2.10) 
The matter-response equation for the perfect-fluid T, (2.5), is 

p = -(p + p)(A+\\) (2.11) 
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3. Vacuum Solution 



With p = p = 0, it is convenient to form the combination \{G* + G x x + ^G y v ), i.e., A + \(A + |A) = 0, 
which is immediately integrated: 

(e A / 2 ), t = c/R 2 (3.1) 
Likewise, the combination ^(G/ + G x x ), i.e., A + A 2 + ^AX = 0, gives 

(3.2) 



R 2 e X/2 

Combining (3.1) and (3.2) leads to a second-order equation for R with solution: 

R 2 = R 2 (t/t )^ (3.3) 

Then 



e 



( C + Cl ) 2 (tA )— (3.4) 



Now, integrating (2.10) and combining the result with (3.3) and (3.4) yields 

i> = % + log(*Ao) (3-5) 

C T" C-y 

In these equations, c, c , and c x are constants of integration which are not all independent. Indeed (2.7) 
and (2.8) (with p = p = 0) are satisfied only if k 2 = (c /c 1 ) 2 = 1 + 4c/c 1 . From this, the final form for the 
vacuum solution becomes: 

e A - Ao = (t/t Q ) 3+^ 



4k 
3Tk 



V> =^o + TT^2 lo g(^o) (3-6) 



Alternatively, a compact form for [} and / is written in terms of the complex parameter e q = (3 + if = R 2 e^: 

e q - q <> = (t/t ) i ^ 1 (3.7) 

As in Einstein's General Relativity, whose corresponding solution (Kasner 1921) can be recovered by setting 
k = and ipQ — mr, n = 0, 2, 4..., the solution is open in the plane of symmetry and closed in the perpendicular 
direction if n 2 < 1. If k 2 > 1, however, the solution is open in all three spatial directions. 

Assuming that the metric is given by g( M „) 5 its signature changes sign every time cost/; = since 
P = R 2 cos ip. This occurs an infinite number of times over a finite interval starting at t = 0, although the 
interval can be made as small as desired by decreasing the NGT parameter k. 

At t — 0, both R and the 4-volume element yj— g go to zero, whereas g n (— e A ) goes to infinity if k 2 < 1 
and to zero if k 2 > 1. The solution would thus appear to contain an essential singularity. This is supported 
by a computer calculation of the Kretschmann scalar, R ol ^ >lv R ol ^ iiv ' , with a symbolic manipulator; the result, 
which we do not reproduce here, contains an additive term that goes like 1/i 4 , along with other terms, also 
in 1/t 4 , but whose numerator oscillates infinitely fast at t = 0. 
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4. Solution with Matter 



In this case, it is best to start from a combination which contains neither p nor p. If we subtract (2.9) 
from (2.8), we obtain after some rearrangement an equation which is also free of ip: 

(log R 2 e~% t + (log ii 2 e A /2) ;t (iog &e-% = (4.1) 

At this point, it is convenient to change the independent variable t to 

f dt , , „ N 

Then we obtain from (4.1) 

(log R 2 e~ X ), vv = (4.3) 

whose solution reads 

R 2 = e *+<x(v-v ) (44) 

whereas (2.10) becomes ip vv = and has for solution: 

tp = a(v — v ) (4.5) 

In (4.4) and (4.5), a and a are arbitrary constants. 

Now we assume an adiabatic equation of state for the matter: p — -fp, with < 7 < 1. From (2.7) and 
(2.8), the combination |(G t * + G x x ), expressed in terms of v, is 

A vv = 8tt(1 - j)R 4 e x p (4.6) 

With all these results, we can write (2.7) in terms of v and of a new dimensionless dependent variable x: 

x v = \/6tt6(1 - 7) {x 2 - rf) ^ (4.7) 

where 

x 2 = R 4 e x p/b + V 2 (4.8) 
and if = {(J 2 + a 2 /3)/32irb, with b an arbitrary constant. The integral of (4.7) is 

x = —rj coth[\/ 6-^677(1 — j)(v — v )] (4.9) 

Now, from the conservation equation (2.11), p{R 2 e x / 2 ) 1+1 is equal to a constant which we identify with 6, 
and this allows us to write (4.8) as x 2 - rj 2 = (i? 2 e A / 2 ) 1-7 . 

Then, with the help of (4.4), we can eliminate exp[3a(l — r y)v/2] between (4.8) and (4.9) to find an 
expression for R as a function of x: 



R = 



Also, by virtue of (4.9), (4.5) becomes 

4 a 1 



4 a 1 t n\ 

lb = — - = -arctanh - (4.11) 
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Finally, from (4.4), there comes 



3(1- T ) 



(4.12) 



To find i as a function of time, we differentiate (4.8) with respect to time and note that v — (R 2 e x ^ 2 ) 1 . 
Then 

V6nb(l -l)(t-t ) = J(x 2 -rj 2 )^ dx (4.13) 

Equations (4.10)-(4.13) constitute the general time-dependent plane-symmetric NGT solution for an 
adiabatic equation of state and — It reduces to the solution obtained in General Relativity (Schiicking 
and Heckmann 1958) in the limit a — > 0. From now on, however, we put equal to zero the constant a 
characterising the anisotropy due to matter in (4.4); any remaining anisotropy is due to NGT. In that case 
we are left with 



R = [X - 7] 



2\ 3(1-7) 



e x = R 2 



4 1 (T]\ 

—= -arctanh — 

V37-1 



p =p (R/R )-^ 



(4.14) 



where now, r\ — a/y/32TTb, a dimensionless constant independent of initial conditions. We can also write the 
first two equations of (4.14) as 



3(1-7) 



(4.15) 



The dependence of x on time is to be found from (4.13). For instance, if 7 = (matter-dominated universe), 
6wb t 2 = x 2 , so that 

R = R (6np t 2 -rf/R 3 ) 1 / s (4.16) 
whereas, if 7 = 1/3 (radiation-dominated universe), R 2 = (x 2 — rj 2 ), with x given by the implicit relation 



V^ot = ^ [x(x 2 v 2 ) 1/2 V 2 log (* + (x 2 - rff' 2 ) 



(4.17) 



It is clear from (4.14) and (4.15) that this solution becomes isotropic in the large x (large time) limit, in 
the sense that / — > and (3 — + e A , so that we recover the Friedmann-Robertson- Walker solution of General 
Relativity It is equally clear that at some time t s < t , there occurs a branch-point at x — rj; R and e A go 
to zero, and the matter density becomes infinite. When 7 = 0, 



8irb 



(4.18) 



One can show from (4.13) that, both for a dust and a radiation universe, the interval between t and t s is 
smaller in the NGT solution than in the corresponding one in GR. 
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5. Early-time Behaviour and Correction to the Age of the Universe 



At t = t Q , the angle ib takes the value 



4 1 
71 7^ 



arctanh-t 



R 3 +V 2 



(5.1) 



whereas ip(t s ) — > — oo. Thus, (3 and / must change signs an infinite number of times between t s and t . The 
zeros of [3 occur when cos tp = 0, at 



£ n = ?7 coth 



^7r(l-7)(2n-l) 



(n = l,2,...) 



(5.2) 



For a dust-filled universe, the first zero of /3 is at = l.7r] ~ a/6b. All other zeros lie between r\ and 1.7r]. 
The time elapsed between t 1 and i is 

.2 \ 1/2 



i ^ — 



1 



- 1.7- 



'■/ 



3/2 



(5.3) 



The behaviour of /? over the interval < t < t x has no obvious physical meaning, but it can be shrunk to 
unobservable size by reducing the value of r\ (or a). Then, if we take t as the present time, expression (5.3) 
can be interpreted as a rough estimate of the age of the universe. 

Now, another way of expressing the age of the universe is with the Hubble parameter, here defined as 
the average rate of expansion of its 3- volume: H — |(log/3e A / 2 ) t , or, if we use the (R, -0) representation, 
H = R/R — |0tan ip. Also, (2.8) can be cast in the form 



R 
R 



8tt 



1 



-p H ib 

3 F \T 



(5.4) 



Insertion of the solution with matter then yields, to first order in (ct/_Rq) 2 : 



12 U 3 



16 



3(l- 7 ) 



8tt / R 



3 \R, 



-3(1-7) 



1 /' " \" / // • 



Rn 



12 \Rl 



1 



16 



3(1-7) 



(5.5) 



From this comes a relation, valid at the present time, between the Hubble parameter, the matter density, 
and (t/Rq-. 

19 / a x - 



Now we can express our "age of the universe" in (5.3) as 



R% 



t t 1 



3H„ 



1 - 



1 a 
15 b 



(5.6) 



(5.7) 



This last form is convenient because both a and b — pR? remain constant throughout the time evolution. 
But no very stringent bound can be inferred from it. For instance, if a/b = 0.1, the length of the matter- 
dominated era is shortened by one percent only. 



In a radiation-dominated universe (7 = 1/3), we have instead 
or, in terms of black-body radiation temperature, 



H 2 = 



1 + 



9 



32™ V T 2 Rl 



T 



(5. 



(5.9) 



where a is the Stefan constant, and a is the same as in the dust solution so as to ensure proper matching of 
the two cases. 



6. Birefringent Universe 

As a nonmetric theory of gravitation, NGT couples matter fields directly to g^ in a way that generically 
violates the Einstein equivalence principle. For instance, the propagation of electromagnetic waves should 
be affected by the coupling between g^ v j and electromagnetic fields. This was first pointed out by Mann 
and Moffat (1981) and detailed calculations were carried out (Will 1989, Gabriel et al 1990, 1991a, 1991b, 
1991c) for gr i0 i 7^ 0, gr^-i = in the static spherically-symmetric case. In the cosmological plane-symmetric 
solution with matter discussed here, we should also find some dependence of the propagation of light on its 
polarisation. 

The coupling of NGT to electromagnetism is written as an electromagnetic Lagrangian density which 
has the general form (Mann et al 1989) 



1 

'16tT 



Z)F av F^ 



(6.1) 



where Y and Z are constants, and T = 1 when g^ = 0. As usual, the electromagnetic field tensor, F^ v , 
can be decomposed into electric components, F i0 = E i7 and magnetic components, F^ — e i j k B k . We also 
take T = sf^g/^J- det g (jUI/) . With (2.1) and R 2 = e x in (4.14), (6.1) can then be written 



1 

8^ 
1 

8^ 



e \\ E \\ + e ± E ± 



E 2 + C-± 

£1 




l)(n-E) 



^)(n-B) 2 



(6.2) 



where 



€ ± = £||/ COS V 

fi ± = jU| cos — 2A sin 2 tp)~^ 



(6.3) 



and n is a unit vector perpendicular to the plane of symmetry. One may think of e as a diagonal tensor 
with nonzero components (e xx , e yy , e zz ) = (ej_, ey , ey), and similarly for a diagonal tensor fi. Our choice of 
T ensures that = fj,^, so that any NGT effect must arise from the (n ■ E) 2 and (n ■ B) 2 terms in C. 
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The Maxwell-NGT equations in momentum space which follow from varying (6.2) are 

k D = 

k B = 

kAE-LoB = 

kAH + LuD = (6.4) 

where we have defined D = e • E, and B = \i ■ H . Note that E and H are not perpendicular to k, but D 
and B arc. To derive a dispersion relation for electromagnetic waves, it is reasonable to treat e and [/, as 
constants over one wavelength. One then obtains 

(e||M||^ 2 - k 2 )B + (l - [k 2 n - k{n ■ k)](n ■ B) + (\ - [B ■ n A k](n A k) = (6.5) 

This leads to two possible speeds for the waves, one for waves with polarisation perpendicular to the 
direction of anisotropy, n, and one for waves polarised in the plane defined by n and k. In the first case, 
n ■ B = 0, and the speed is 

c_l = l(l-(l-cosi/0sin 2 6>) 1/2 (6.6) 

where 9 is the angle between the direction of anisotropy and the direction of propagation of the waves. 

On the other hand, when their magnetic field lies in the plane that contains n and k, waves propagate 
at speed 

C|| = - (l - [(1 - cos V) + (2X - 1) sin 2 VV cos ip] sin 2 9 J (6.7) 

There are two contributions to (6.6) and (6.7) which should be carefully distinguished. One arises simply 
from the fact that the waves propagate on a NGT plane-symmetric background. Indeed, consider a source 
situated at a distance r from the Earth, with the vector r making an angle 9 with n (or x). From (2.1), the 
equation for light propagation is 

dt 2 - e x dx 2 - [3{dy 2 + dz 2 ) = (6.8) 

Then, in the (R,ip) representation, and using e A = R 2 from (4.14), the speed of the light emitted by the 
source as measured by an observer on Earth would be 

% = - [1 - (1 - cos iP) sin 2 9] 1/2 (6.9) 
dt R 1 J 

which is precisely Cj_. In fact, only the term proportional to sin 2 tp in is a consequence of the direct 

coupling between NGT and electromagnetism. 

The resulting effect bears a lot of similarity to birefringence in crystals with one optical axis. The 

polarisation vector of an electromagnetic wave rotates along its trajectory in a way that depends on the 

difference between the two speeds. Assuming that ip is very small all along the trajectory in a dust-filled 

universe (7 = 0), this is: 

Cm — Cj_ w \;{X — i) sin 2 tp sin 2 9 

" R 

"5<*-i>(£)"(KSF? <«»» 

The only parameter-free prediction here is that waves propagating along the anisotropic direction will expe- 
rience no rotation of their polarisation vector. 
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7. Conclusion 

Equations (3.6) give the solution to the vacuum homogeneous, plane-symmetric case in NGT, whereas 
equations (4.14) represent the equivalent solution in the presence of a perfect fluid. Both exhibit zeros in 
g 22 and in gr 23 ]- The actual physical meaning of this behaviour is not known at present. In the solution with 
matter, the NGT parameter a can be chosen so that the zeros occur early enough in the evolution of the 
universe for them to be essentially indistinguishable from the initial singularity. Thereafter this anisotropic 
universe approaches a FRW universe in such a way that currently well-established features of its early history 
(such as helium production) do not change significantly. 

It is interesting that these signature reversals do not occur in a complex Hermitian version of NGT, 
with g^ v = g( M „) + e S[ f _ tu ]' w here e 2 = —1. Taking [3 = R 2 cosh?/' and / = R 2 sinh^>, the solution (3.6) now 
becomes 

-2(1 + k 2 ) 

e A ~ A ° = (*/*o) 3 ~" 2 
R/R = (t/t )^ 

V> = ^o + — ^2 lo g(*Ao) (7-1) 

4(1 + k) -8k 

with k a real constant. Therefore, (3 goes like t 3 -*' 2 (1 + t 3 -* 2 ) and does not go through zero. If k 2 > 1, it 
even has a minimum at t = (-^tt) 8k . 

In the case with matter, the solution (4.14) becomes 

R = (x 2 + V 2 )^ 

4 1 /77\ 

w = —= -arctan — 

V37-1 

e x = R 2 

p =p a (R/R a )- 3 ^ (7.2) 

where r\ = a/^/32nb is real, and x is still obtained from (4.13). As t — > (or x — > 0), tp — > (2tv/V3)/ (7 — 1), 
and (3 — > ^ 4 / 3 ( 1 ~t) cosh[(27r/v / 3)/(l — 7)], again without any oscillation through zero. 

It has been asserted (Kelly and Mann 1986, Mann and Moffat 1982) that only the hypercomplex structure 
leads to a linearised NGT that is free of ghost poles. Yet, as mentioned in the introduction, hypercomplex 
NGT does have ghost poles generically if the modes are not somehow decoupled to all orders. Thus, 
it may be premature to reject the Hermitian structure, especially since, at least in the cases discussed here, 
the solutions do not exhibit the puzzling signature reversals of the hypercomplex theory. Whether their 
singularity structure differs from that encountered in GR, as suggested by (7.1) and (7.2), can only be 
established by a more thorough investigation. 

Finally, an intriguing possibility is that some small residual anisotropy may affect the propagation of 
light in such a way that its polarisation vector rotates as it travels over cosmological distances. Such an effect 
has been reported recently (Nodland and Ralston 1997a, 1997b, 1997c), but the validity of these observations 
has been disputed (Carroll and Field 1997, Eisenstein and Dunn 1997, Leahy 1997, Loredo et al 1997, Wardlc 
et al 1997). If they were confirmed, NGT would provide a possible explanation; if not, an upper bound could 
be obtained on the size of the coupling of NGT to electromagnetism. 
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